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Analysis  of  the  absorption  spectrum  of  rnby  at  high  pressures 


Surinder  M.  Sharma4  and  Y.M.  Gupta 
Shock  Dynamics  Laboratory 
Department  of  Physics 
Washington  State  University 
Pullman,  WA  99164-2814 


It  is  shown  that  at  high  pressures,  results  of  absorption  measurements  are  inconsistent  with 
R-line  fluorescence  data.  The  observed  constancy  of  Rj-Rj  splitting  under  hydrostatic  conditions 
is  supported  by  crystal  field  analysis.  The  present  findings  suggest  the  need  for  absorption  meas¬ 
urements  in  ruby  crystals  under  better  defined  stress  conditions. 


1988  PACS  numbers:  62.50.  +p,  78.20.Dj,  78.55.Hx 
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Ruby  spectra  have  been  the  subject  of  numerous  experimental  and  theoretical  investigations 
due  to  the  utility  of  ruby  as  a  laser  material  and  as  a  high  pressure  sensor  in  diamond  anvil  cell 
studies.1-7  At  ambient  conditions,  various  details  of  the  spectrum  consisting  of  multiplets  of  R, 
R'  and  B  lines,  and  U  and  Y  absorption  bands  can  be  quantitatively  understood  with  symmetry 
based  crystal  field  theory  approach.5,8-10  Significant  crystal  field  parameters  for  description  of 
ruby  spectrum  turn  out  to  be  the  following:5  Racah  parameters  B  and  C,  representing  electron- 
electron  repulsion  of  d-electrons  of  Cr3+;  10  Dq,  characterizing  the  strength  of  crystal  field  of 
octahedral  symmetry;  V lrigonai,  the  crystal  field  of  only  C3  symmetry  which  is  the  local  site  sym¬ 
metry  around  Cr3+  in  the  corrundum  lattice.11-13  The  trigonal  field  is  characterized  by  the  magni¬ 
tude  of  two  matrix  elements  of  this  field  (between  the  bases  of  symmetry  adapted  sets  {t2s}  and 
(ej)  namely  v  =  -3<t2„x+  |  V trig  |  t2„x+>  and  v  =  <t2jx+|  vtrit  I  es«  +  >.  as  defined  by  Pryce 
and  Runciman.14  Fine  multiplet  structure  of  energy  levels  emerges  from  the  simultaneous  presence 
of  trigonal  field  and  spin-orbit  interaction  and  can  be  explained  quantitatively.5,8  9 

Attempts  to  understand  variations  in  ruby  spectrum  due  to  hydrostatic  pressure,  in  terms  of 
crystal  field  approach,  have  aimed  at  broad  correlations.15,16  The  fine  multiplet  structure  of  lines 
and  bands  is  completely  ignored  —  thereby  approximating  the  crystal  field  to  have  octahedral 
symmetry  only.  Variations  in  crystal  field  parameters,  in  the  assumed  octahedral  symmetry,  are 
represented  by  two  independent  ‘scale  parameters’,15  16  which  are  determined  by  the  best  fit  to 
appropriately  averaged  experimental  data.  Such  a  procedure,  by  its  very  nature,  overemphasizes 
the  consistency  of  pressure  dependent  changes  in  the  spectrum.  For  example,  over  80  percent  of 
the  measured  shift  of  the  Y  band  (representing  47’1— *4.42  transition)  is  due  to  changes  in  10  Dq 
and,  therefore,  a  considerable  variation  in  the  pressure  dependence  of  Racah  parameters  with  pres¬ 
sure  will  not  impair  the  quality  of  agreement. 

In  this  communication  we  analyze  the  consistency  between  the  pressure  variation  of  the 
absorption  bands,17  and  Ri,  R2  lines18  in  detail  i.e.  beyond  the  octahedral  approximation  in  the 
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crystal  field  analysis.  The  consistency  analysis  is  restricted  to  this  part  of  the  spectrum  for  two 
reasons.  First,  no  other  lines  have  such  a  well  defined  spectrum  as  R-lines;19  second,  this  restricted 
analysis  suggests  the  need  for  experimental  measurements  of  the  absorption  bands  under  well 
characterized  conditions  of  stress,  and  with  improved  statistical  precision.  An  accurate  determina¬ 
tion  of  the  details  of  the  ruby  spectrum  at  high  pressure  would  permit  a  detailed  crystal  field 
analysis  similar  to  that  performed  by  Macfarlane9  at  ambient  conditions. 

The  results  describing  the  pressure  variation  of  the  absorption  spectrum,  which  enter  as  input 
to  our  calculation,  have  been  taken  from  the  polarized  light  absorption  measurements  of  Stephens 
and  Drickamer.17  They  measured  the  absorption  bands  up  to  120  kbar  and  used  a  crystal  field 
description  to  derive  10  Dq,  B  and  trigonal  field  as  a  function  of  pressure.  These  authors  assumed 
that  both  the  bands  viz  U  (4T2)  and  Y  (^Ti)  split  by  the  same  amount  v/2.  Subsequent  analysis 

of  Macfarlane9  showed  that  the  splitting  of  Y  band  is  given  by  ~  +  v.  Hence,  only  the  experi- 

mental  results  of  U  (4T2)  band  should  be  used  to  estimate  changes  in  trigonal  field.  Further,  by 
comparing  the  splitting  of  U  and  Y  bands,  we  conclude,  that  variation  in  v  with  pressure,  within 
experimental  errors,  is  zero.  The  observed  variations  of  absorption  bands  from  Ref.  17  are  sum¬ 
marized  as  follows,  (i)  Racah  parameter,  B,  decreases  as  .205  cm- ‘/kbar  (»)  io  Dq  increases  at  a 
rate,  10.9  cm'’/ltbar  (iii)  trigonal  field  parameter  v  shows  little  or  no  change  up  to  «75  kbar, 
beyond  which  it  increases  rapidly.  We  now  use  these  results  to  calculate  the  variations  in  R„  R, 
lines  with  pressure.  The  calculation  of  the  red  shift  is  followed  by  the  calculation  for  R-lines  split¬ 
ting. 

Reduction  in  Racah  parameter  B  with  pressure  implies  a  red  shift  of  the  R  lines  and  is,  there¬ 
fore,  in  qualitative  agreement  with  experimental  observations.18  If  there  is  no  change  in  the  trigo¬ 
nal  field  with  pressure,  then  the  rate  of  the  red  shift  with  pressure  can  be  deduced  by  diagonalizing 
the  4x4  matrix  representing  Coulomb  interaction.5  With  the  above  mentioned  observed  values,17 
the  red  shift  of  R  lines  is  given  by  4.26  cm-1/kbar  (Assume  C/B  =  5.1245  to  get  position  of  R  lines 
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at  ambient  conditions].  This  is  to  be  compared  with  the  measured18  red  shift  of  .76025  cm-1/kbar; 
this  is  smaller  by  a  factor  of  five  when  compared  to  the  calculated  shift.  A  calculation  based  on 
sensitivity  analysis,  after  Macfarlane,9  which  includes  configuration  interaction,  confirms  that  the 
calculated  red  shift  is  about  a  factor  of  five  larger  than  the  observed  value.  This  implies  that  the 
variation  of  Racah  parameter,  obtained  from  absorption  experiments,  is  in  significant  disagreement 
with  the  red  shift  of  R  lines. 

As  pointed  out  earlier,  coexistence  of  the  trigonal  field  and  spin-orbit  interaction  bring  about  split¬ 
ting  of  2E  giving  rise  to  the,  Rt,  R2  lines.  Within  second  order  perturbation  theory,  the  splitting 
between  Ri ,  R2  lines  can  be  written  as 

4 

3Vf 

a(Ri  _  r2)  =  - 4 -  •{!) 

where  f  represents  the  strength  of  the  spin-orbit  interaction.  Eq.  (1)  can  be  used  to  predict  the 
change  in  splitting  with  changes  in  the  trigonal  field.  We  shall  assume  that  spin-orbit  interaction 
does  not  change  with  pressure  —  an  approximation  that  is  consistent  with  an  analysis  of  the  life¬ 
time  variation  of  Rt  line  with  pressure.20  Alternatively,  we  can  use  the  following  result  due  to 
Macfarlane21  to  evaluate  the  changes  in  R-line  splitting  due  to  variations  in  the  trigonal  field 

A(RX  -  R2)  =  D0  -  0.027 -  0.008^ '  (2) 

Here,  D„  represents  splitting  at  P=0,  and  vp  and  vp>  are  pressure  dependent  changes  in  v  and 
For  the  same  input,  Eq.  (1)  gives  a  larger  splitting  than  Eq.  (2)  and  so  we  term  the  results  of 
Eq.  (l)  as  upper  and  that  of  Eq.  (2)  as  lower  bound.  In  principle  Eq.  (2)  should  be  more  accurate, 
as  it  includes  configuration  interaction.  We  have  used  both  equations  because  experimental 
results19  of  relative  movements  of  R  and  R'  lines  do  not  confirm  the  predictions  of  Macfarlane’s 
calculation.9 
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We  now  calculate  the  R-line  splitting  due  to  the  trigonal  field  changes,  reported  on  the  basis 

of  the  absorption  measurements.17  We  take  €(2 T2)  -  €(2£)  from  Ref.  9  and  f  =  186.1  cm'1,  to 
reproduce  R-line  splitting  at  ambient  pressure.  Results  for  a  few  representative  pressures  are  given 
in  Table  I. 

Therefore,  the  absorption  data  of  Stephens  and  Drickamer17  imply  no  change  in  splitting  up 
to  about  75  kbar  and  then  a  large  change.  At  120  kbar,  the  splitting  increases  by  approximately 
30%  compared  to  that  at  75  kbar.  However,  the  experimental  results  for  R-lines,  consisting  of  38 
sets  of  data  up  to  150  kbar,  tested  for  statistical  consistency,18  are  incompatible  with  the  results 
above  75  kbar.  The  experimentally  observed  rate  of  change  of  splitting  is18  .0025  ±  0.004 
cm-1/kbar.  This  result  implies  that  within  the  experimental  accuracy,  there  is  no  change  in  split¬ 
ting  up  to  150  kbar,  in  contrast  to  the  results  shown  in  Table  I. 

Experimental  results  of  R-lines,  aside  from  better  precision  and  mutual  statistical  compatibil¬ 
ity,18  are  also  consistent  with  crystal  field  analysis  as  shown  below.  We  use  the  point  ion  approxi¬ 
mation  to  estimate  changes  in  the  crystal  field.  Although  this  approximation  is  known  to  be  inade¬ 
quate  for  calculation  of  crystal  field  parameters  at  ambient  conditions,10'22  it  gives  a  reasonable 
estimate  of  stress  induced  changes.21,23,24. 

As  the  site  symmetry  of  Cr3+  site  is  C3,  the  potential  function  for  an  electron  can  be  written3 


V=  flg  +  B^rCoaO  +  B%r2(ZCoa28  -  1)  +  5§r3(5 Coa39  -  3Cos9) 

+B\r3Sins9Co8Z<f>  +  Clr3Sin39Sin34> 

+Boiri(35Cos*0  -  30 Coa29  +  3)  +  B3iriSin39Cos9Cos34>  +  ■  ■  ■  (3) 

where  5JJ1  and  CJJ*  are  the  lattice  sums  over  real  and  imaginary  part  of  the  functions  of  the  type 


R 


»+i 


Ylm  as  defined  and  calculated  by  McClure3  at  ambient  conditions  [the  constant  j  depends 


on  the  position  of  the  atom  with  respect  to  the  symmetry  axis  and  R  represents  lattice  vectors). 
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The  trigonal  field  parameters  v  and  v  can  now  be  represented  in  terms  of  B™  as  follows3  23 


6  rjO  2  160 

.  =  7Jfi<r«>  +  — 


5?  + 


20y/2 


<r*> 


-2V2  „n  .  a.  .  40 V2" 

V  =  — - — 5j<r  >  + 


52  + 


*2 


20%/2 


<r4> 


(4) 

(5) 


7  63 

We  assume  that  hydrostatic  pressure  results  in  an  isotropic  compression  of  the  rhombohedral 
cell25-28  of  ruby  and  that  all  interatomic  distances  change  uniformly  and  proportionally  to  unit  cell 
dimensions  a.  This  assumption  seems  to  be  reasonable  up  to  as  3%  compression  («  100  k bar). 29 

SB Sa 

Then  it  can  be  shown  that  —  =  (n-t-1) —  for  small  Sa/ a.  Such  an  approximation  gives  rea¬ 
sonable  agreement  with  experiments.30  The  changes  in  v  and  v  under  pressure  cam  then  be  writ¬ 
ten  as  follows 


Sv  —  — 


c  §. a 

0V  =  - 


Sa 


18  *^0  2  800  4 

—  5^<rz>  +  <r4> 

7  63 


52 


6y/2 


200 \/2 


<r2>B$  -  5$  + 


52  ' 
20\/2 

g?  11 

20^2 
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Taking<r4>/<r  2>  =  3.0732 

a3  estimated  by  Schawlow  et  al.  24  from  Hartree-Fock  solutions  for  a  free  Cr3+  ion,  we  get 

Sa 


(6) 

(7) 

(8) 


Sv  =  —<r2>.  101803 
a 


(9) 

Sv  =  — <  r2  >  .256326  (10) 

a 

Eqs.  (9)  and  (10)  give  a  negative  Sv  and  a  positive  Sv  consistent  with  earlier  conclusions.21 
Substituting  these  values  in  Eq.  (2)  gives 


A(R!  -  R2)  =  D0  +  6.981  x  10~4<r2>  — 


Sa 


For  ruby  —  =  -1.28  x  10~45,  where  P  is  pressure  in  kbar.25  27  This  gives  a  splitting  of 

a 


A(Rt  -  R2)  =  D0  -CXP 


(12) 


Where  C\—  8.935  x  10-8<r2>  ev/kbar.  The  value  of  <r2>  has  been  estimated  to  be  .4048 
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A 2  with  Hartree-Fock  wavefunction  and  1.38  A 2  with  Slater’s  wavefunction  by  Schawlow  et  al.24 
These  give  [Ci\aF  =  2.92xl0-4  cm_1/kbar  and  [C,1]Sj(1<e,  =  9.95x  10~4  cm'V^bar  implying  that 
the  splitting  changes  at  the  rate  .001  cm-1/kbar  or  less  and  is  negligible.  This  finding  is  consistent 
with  observed  results.  At  compressions  larger  than  3%,  local  compressibility  around  Cr3+  may  be 
smaller  than  the  bulk  compressibility29  of  AljOj.  Then  the  above  analysis  would  give  an  upper 
limit  on  variation  of  trigonal  field  within  the  point  ion  approximation. 

On  the  basis  of  the  above  analysis,  we  conclude  that  the  R-line  data  are  consistent  with  crys¬ 
tal  field  model.  The  high  pressure,  absorption  results  reported  in  Ref.  17  are  not  consistent  with  R 
line  data,  on  both  counts  —  the  red  shift  as  well  as  splitting  of  Rj,R2  lines;  the  splitting  results  are 
consistent  to  75  kbar.  We  conjecture  that  the  inconsistency  between  the  absorption  and  the  R-line 
results  may  arise,  in  part,  due  to  the  possibility  of  non-hydrostatic  stresses.  It  is  known  that  non¬ 
hydrostatic  stresses  due  to  either  uniaxial  stress51  or  uniaxial  strain32  bring  about  changes  in  the 
trigonal  field.  Hence,  the  existence  of  non- hydrostatic  stresses  beyond  75  kbar  in  the  experimental 
configuration  used  in  Ref.  17  may  explain  the  increase  in  the  trigonal  field.  Experimental  data  are 
needed  to  verify  this  conjecture.  In  addition,  to  accurately  determine  the  adequacy  of  the  crystal 
field  description  for  modeling  ruby  spectrum  at  high  pressures,  there  is  a  need  for  high  precision 
absorption  data  under  well  defined  hydrostatic  pressures. 

Before  concluding,  we  would  like  to  make  another  observation.  Ma  et  al.16  assume  an  effec¬ 
tive  interaction  between  2E  and  2T^  states  which  is  an  increasing  function  of  trigonal  field.  To 
explain  the  experimental  observations  of  the  relative  shift  of  R  and  R'  lines  they  require  a  large 
variation  of  the  trigonal  field  with  pressure.  This  is  in  contrast  with  our  conclusion  that  the  trigo¬ 
nal  field  does  not  change  significantly  with  hydrostatic  pressure.  Therefore,  the  assumed  interac¬ 
tion  of  Ma  et  al.16  requires  a  reconsideration.  In  fact,  V,rig  does  not  couple  the  2 E  and  2T i  states 
to  the  first  order .1,s  Significant  configuration  interaction  is  through  2  T2  and  is  therefore  of  second 
order  in  perturbation.  Ma  et  al.16  require  a  fairly  strong  interaction  between  2E  and  2  T t  states 
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and,  further,  it  needs  to  increase  rapidly  with  pressure.  We  feel  that  such  an  effective  interaction 
actually  emerges  because  of  increased  mixing  of  oxygen  p  orbitals  with  that  of  Cri+  d  or bi tails. 13 
Such  a  mixing  is  expected  to  be  mainly  responsible  for  the  removal  of  accidental  degeneracy  of  2  E 
and  27\  levels  which  exists  within  first  order  perturbative  treatment  with  d  orbitals.5  The  sarnie 
could  well  be  the  cause  of  disagreement  of  observed  sensitivity  of  R'  lines19  with  crystail  field 
theoretical  predictions.9 

Discussions  with  P.  Horn  and  J.  Burt  about  the  ruby  spectra  have  been  helpful  to  us.  This 
work  was  supported  by  ONR  Contract  N0O-014-86-K-0307. 
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TABLE  I: 


Variation  in  R1-R2  splitting  due  to  observed17  changes  of  the  trigonal  field.  Eq.  (i) 
gives  the  upper  and  Eq.  (2)  gives  the  lower  bound. 


Pressure  (kbar) 

Upper  bound  (cm  *)  1 

Lower  bound  (cm1) 

0 

29.6 

29.6 

1 

75 

29.6 

29.6 

92 

33.4 

32.4 

102 

36.4 

34.6 

120 

39.7 

37.0 

